Magnetic switching of a single molecular magnet due to spin-polarized current 
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Magnetic switching of a single molecular magnet (SMM) due to spin-polarized current flowing 
between ferromagnetic metallic leads (electrodes) is investigated theoretically. Magnetic moments 
of the leads are assumed to be collinear and parallel to magnetic easy axis of the molecule. Elec- 
trons tunneling through the barrier between magnetic leads are coupled to the SMM via exchange 
interaction. The current flowing through the system as well as the spin relaxation times of the SMM 
are calculated from the Fermi golden rule. It is shown that spin of the SMM can be reversed by 
applying a certain voltage between the two magnetic electrodes. Moreover, the switching may be 
visible in the corresponding current-voltage characteristics. 

PACS numbers: 75.47.Pq, 75.60.Jk, 71.70.Gm, 75.50.Xx 



I. INTRODUCTION 

Although first synthesized in the 1980s, single molec- 
ular magnets (SMMs) 1 did not get much attention un- 
til the beginning of the 1990s, when their unusual mag- 
netic properties were discovered 2 . Owing to large spin 
and high anisotropy barrier, SMMs in a time dependent 
magnetic field were shown to exhibit magnetic hystere- 
sis loops with characteristic steps caused by the effect 
of quantum tunneling of magnetization. Current in- 
terest in SMMs is a consequence of recent progress in 
nanotechnology, which enables to attach electrodes to 
a single molecule and investigate its transport proper- 
ties 3 ~ 7 . Physical properties of SMMs and their nanoscale 
size make them a promising candidate for future applica- 
tions in information storage and information processing, 
as well as in various spintronics devices 8 . 

Magnetic switching of a SMM due to quantum tunnel- 
ing of magnetization in a magnetic field varying linearly 
in time was considered theoretically long time ago and 
was also studied experimentally 2 . From both practical 
and fundamental reasons it would be, however, interest- 
ing to have a possibility of switching the SMM without 
external magnetic field. Such a possibility is offered by 
a spin polarized current. As it is already well known, 
spin polarized current can switch magnetic layers in spin 
valve structures, like for instance magnetic nanopillars 9 . 
The main objective of this paper is just to investigate 
theoretically the mechanism of SMM's spin reversal due 
to spin polarized current. 

As a simplest system for current-induced molecular 
switching we consider a SMM embedded in the barrier 
between ferromagnetic electrodes (called also leads in the 
following). When voltage is applied, the charge current 
flowing in the system is associated with a spin current. 
In this paper we show that this spin current can lead to 
magnetic switching of the SMM, when the voltage sur- 
passes a certain threshold value. Moreover, when bias in- 
creases (linearly) in time, the switching can be observed 
in the corresponding current- voltage characteristics as an 
additional feature (dip or peak) in the current. 



It is worth to note that spin polarized transport 
through artificial quantum dots attached to ferromag- 
netic leads was extensively studied in recent few years, 
mostly theoretically 10-12 , though some experimental 
data are already available 13 . However, investigations of 
spin polarized electronic transport through molecules, 
and particularly through magnetic ones, are in early 
stage of development. 

The paper is organized as follows. In section 2 we 
present the model Hamiltonian assumed to describe a 
molecule interacting with magnetic leads. Theoretical 
analysis of electric current flowing through the system 
under consideration is carried out in section 3. Numeri- 
cal results on electric current and magnetic state of the 
molecule are presented and discussed in section 4. 



II. MODEL 

We consider a model magnetic tunnel junction which 
consists of two ferromagnetic leads separated by a non- 
magnetic barrier, with a SMM embedded in the barrier. 
Electronic transport in the system occurs owing to tun- 
neling processes between the leads. However, the tunnel- 
ing electrons can interact with the SMM via exchange in- 
teraction, leading to spin switching of the molecule. For 
simplicity, we will consider only collinear (parallel and 
antiparallel) configurations of the leads' magnetic mo- 
ments. In addition, magnetic moments of the leads are 
parallel to the magnetic easy axis of the SMM, as shown 
schematically in Fig. 1(a). 

For the sake of simplicity we assume that the spin num- 
ber S of the molecule is constant, i.e. it does not change 
when current flows through the system. This also means 
that the charge state of SMM is fixed and only pro- 
jection of the molecule's spin on the quantization axis 
(anisotropy axis) can be changed due to the current. In 
addition, we restrict the following discussion to the case 
of weak coupling between the molecule and electrodes. 

The full Hamiltonian of the system under consideration 
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FIG. 1: (color online) (a) Schematic picture of the system un- 
der consideration for two collinear configurations of the elec- 
trodes' spin moments; parallel (blue arrows) and antiparallel 
(red arrows). Dashed lines represent the two possible tunnel- 
ing processes: direct tunneling (the top line) and tunneling 
with scattering on the SMM's spin due to exchange interac- 
tion (the bottom line), (b) Energy levels corresponding to 
different spin states of the SMM. The grey dot represents the 
initial spin state, | — S), of the molecule. 



takes the form 



H = 7~Lsmm + 7~Ll + 7~Lr + Ht- 



(1) 



The first term describes the SMM and is assumed in the 
form, 



H 



SMM 



= -DSi 



(2) 



where S z is the z component of the spin operator, and 
D is the uniaxial anisotropy constant. Although Eq. (2) 
represents the simplest Hamiltonian of a free SMM, it is 
sufficient for the effects to be described here. The next 
two terms describe ferromagnetic electrodes, 



ka 



(3) 



for q = L (left lead) and q — R (right lead). The elec- 
trodes are characterized by conduction bands with the 
energy dispersion e£ , where k denotes a wave vector 
and a is the electron spin index. In Eq.(2) and aft 
are the relevant annihilation and creation operators, re- 
spectively. 

The last term of the Hamiltonian Ti stands for the 
tunneling processes 14-16 , 
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The first term in the above equation describes tunnel- 
ing with simultaneous interaction of tunneling electrons 
with the SMM via exchange coupling, with J q q ' be- 
ing the relevant exchange parameter. In a general case 
Jl,l ^ Jr,r Jl.r = Jr,l- In the following, however, 
we assume symmetrical situation, where Jl,l = Jr,r = 



Jl,r — Jr.l = J- The second term of Eq. (4) describes 
direct tunneling between the leads, with Td denoting the 
corresponding tunneling parameter. Apart from this, S 
is the SMM's spin operator, and cr — (a x , a v , a z ) is the 
Pauli spin operator for conduction electrons. We assume 
that both, Td and J, are independent of energy and po- 
larization of the leads. Additionally, Td and J are nor- 
malized in such a way that they are independent of the 
size of electrodes, with N q (q = L, R) denoting the num- 
ber of elementary cells in the q-th electrode. 

The electric current flowing in the system is deter- 
mined from the Fermi golden rule, 16 , 

mm' kk'a/3 I, 



(L'kam <->■ Rk' f3m! 



(5) 



where e is the electron charge (for simplicity we assume 
e > 0, so current is positive when electrons flow from 
left to right), /(e) is the Fermi-Dirac distribution, P m is 
the probability to find the SMM in the spin state |m), 
and Wftk"p£i> is the ra te of electron transitions from the 
initial state {L'kam} to the final one {Rk 1 f3m'}. 

Up to the leading terms with respect to the coupling 
constants Td and J, the current is given by the formula 



I = 



DfDf 



DfD?)V 



^e 2 [|T d | 2 + |J| 2 <S 2 )]( 



DfDfA- v (m)dD(-r)2m + 1) + r/eV^j 
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(6) 



Here, D% is the density of states (DOS) at the Fermi 
level in the g-th electrode for spin cr, (S* 2 ) = Y^m m2 Pm, 
and V is the voltage between the leads, eV = [i^ — Mfl> 
with hl and hr denoting the electrochemical potentials 
of the leads. Finally, A±(m) = S(S+ 1) — m(m± 1), and 
C(e) = e[l - exp(-e/3)] ^ with = k B T. 



III. THEORETICAL ANALYSIS 

To calculate electric current from Eq. (6) we need to 
know the probabilities P m . To find them, we assume the 
initial state of the SMM's spin to be | — S) , as indicated 
in Fig. 1(b). By applying a sufficiently large voltage, one 
can switch the molecule to the final state \S). The rever- 
sal process takes place via the consecutive intermediate 
states; | — S+ 1), . . . , \ S — 1). In the following we assume 
that the voltage applied to the system grows linearly in 
time, V = ct, where c denotes the velocity at which the 



voltage is increased. This allows to observe switching di- 
rectly in the current flowing through the system when the 
voltage surpasses a critical value. The probabilities P m 
can be then found from the following master equations: 
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+ lm+l P m+l + 7m-l P ' 



(7) 



m — 1 5 



. cP-s = ~7-s p -S + 7-s+i p -s+i: 



for -S < m < S and P defined as P = dP/dV. The 
transition rates ^jm are given by 

^ = ^\J\ 2 A ± (m) 



\p^Df c(-D(±2m + 1) ± eV 
+ Df L>f ((D(±2m + 1) T eV 



DfDf 



D?Df 



C(D(±2m + l))}. (8) 



The relevant boundary conditions are: P_,s(V = 0) = 1 
and P m (V = 0) = 0, for m 7^ -5. 



IV. NUMERICAL RESULTS AND DISCUSSION 

Numerical calculations have been performed for an oc- 
tanuclcar iron(III) oxo-hydroxo cluster of the formula 
[Fe 8 2 (OH) 12 (tacn) 6 ] 8+ (shortly Fe 8 ), whose total spin 
number is S = 10. The anizotropy constant is D — 0.292 
K 17 , and we assume that J w w 100 meV. Further- 
more, both the leads are assumed to be made of the same 
metallic material, with the elementary cells occupied by 
2 atoms contributing 2 electrons each. The density of 
free electrons is assumed to be n w 10 29 m~ 3 . The elec- 
trodes are characterized by the polarization parameter 
P" = {D\ - Dl)/{D q + + D g _), where D g +( _ } denotes 
the DOS of majority (minority) electrons in the q-th elec- 
trode. The temperature of the system is assumed to be 
T = 0.01 K, which is below the blocking temperature 
T B = 0.36 K of Fe 8 . 

Let us begin with the case where both electrodes are 
nonmagnetic. In Fig. 2(a) we show the average value of 
the z component of the SMM's spin, (S z ) = S n =— S n -Pm 
and the charge current I. The spin reversal is not 
found, though the current affects the SMM's spin for 
V exceeding the threshold voltage determined by the 
anisotropy constant D (energy level separation). At this 
voltage, transport associated with spin-flip of the con- 
duction electrons becomes energetically allowed, exciting 
the molecule to the spin state | — S + 1). As the voltage 
is increased further, the different SMM's spin states \m) 
become equally probable and (S z ) — > 0. 

The situation becomes significantly different when the 
electrodes are ferromagnetic, and tunneling processes are 
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FIG. 2: (color online) (a) The mean value of the SMM's spin 
(S z ) and the current I flowing through the system as a func- 
tion of voltage V in the case of nonmagnetic electrodes, and 
for c = 10 kV/s and T = 0.01 K. The inset shows schemati- 
cally the density of states and the allowed tunneling processes 
(the solid arrows correspond to electrons tunneling without 
spin reversal, whereas the dashed arrows correspond to elec- 
trons tunneling with simultaneous spin-flip), (b) Schematic 
representation of the DOS and tunneling processes for both 
parallel and antiparallel magnetic configurations, shown in 
the case when the left electrode is made of a half-metallic fer- 
romagnet while the right one is either a typical 3d ferromagnet 
or a half- metallic ferromagnet. 



strongly spin dependent. The following discussion is lim- 
ited to the most interesting situation, when one (say the 
left) electrode is a half-metallic ferromagnet with fully 
spin-polarized electrons at the Fermi level, P L = 1. The 
second electrode can be either nonmagnetic, or typical 3d 
ferromagnet, or even half- metallic. Tunneling processes 
with and without spin flip are indicated schematically 
in Fig. 2(b). The corresponding transport characteristics 
and the average value of S z are shown in Fig. 3 for both 
parallel and antiparallel magnetic configurations of the 
leads, and for various spin polarizations of the right elec- 
trode. The complete reversal of the SMM's spin becomes 
now possible, independently of the magnetic polariza- 
tion of the right electrode. Starting with the spin state 
I — S) at zero bias, one arrives at the state \S) when the 
bias voltage surpasses the threshold value. Moreover, the 
switching leads to some features in the tunneling current. 
We note that switching also takes place for < P L < 1, 
but the switching time becomes longer. 

In the parallel configuration, Figs 3(a)-(c), the rever- 
sal process can be observed as a dip in the current, which 
becomes more pronounced when P R — > 1. The dip cor- 
responds to the voltage range where the SMM's spin re- 
versal process takes place. It begins at the same voltage, 
V ~ 0.48 mV, which corresponds to the energy gap be- 
tween the SMM's spin states | — S) and | — 5+1) (approx- 
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FIG. 3: (color online) The mean value of the SMM's spin (S z ) 
and the current / flowing through the system as a function of 
voltage V for indicated polarization parameters in the parallel 
(P) (a)-(c) and antiparallel (AP) (d)-(f) configurations, and 
for c = 10 kV/s and T = 0.01 K. Part (g) corresponds to the 
case with one electrode being nonmagnetic. 
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FIG. 4: The mean value of the SMM's spin (S z ) and the 
current J as a function of the voltage V for various voltage 
sweeping speeds c in the parallel magnetic configuration of the 
system. The numerical results are for T = 0.01 K, P L — 1 
and P R = 0.9. 



imately 5.55 K in the case considered). Because the en- 
ergy gaps between the higher spin states are smaller, this 
energy is the activation energy for the current induced 
switching. Below the threshold voltage only direct tun- 
neling (described by the second term in Eq.(4)) and the 
non-spin-flip part of the first term in Eq.(4) contribute 
to charge current. When the voltage activating spin re- 
versal is reached, some of the tunneling electrons can flip 
their spins due to exchange interaction with the molecule, 
and this leads to spin reversal of the SMM. As a result 



\(S Z )\ becomes reduced. This leads to partial suppres- 
sion of the non-spin-flip contribution to current from the 
first term in Eq.(4). Instead of this, a spin-flip contribu- 
tion becomes nonzero. However, the latter contribution 
is small as it involves DOS in the minority electron band 
of the right electrode, and cannot compensate the loss 
of current due to the non-spin-flip tunneling (which in- 
volves DOS for majority electrons). This leads effectively 
to a dip in the current, which occurs in the voltage range 
where spin switching of the SMM takes place. The dip 
disappears when spin of the SMM is completely reversed. 
The broadening of the dip, in turn, stems from the fact 
that as P R — > 1, the transition times, l/jm^ + \ become 
longer and longer (see Eq. (8)), and the time required for 
complete SMM's spin reversal becomes longer as well. 
This also makes the dip more pronounced. 

The situation is significantly different in the antipar- 
allel configuration, Figs 3(d)-(f). Instead of the dip in 
current, there is now a peak in the voltage range where 
the switching takes place. This is because now the role 
of spin minority and spin majority electron bands in the 
right lead is interchanged. Additionally, the current flow- 
ing through the system tends to when P R — ► 1 (perfect 
spin valve effect), except for a small voltage range where 
the reversal of the SMM's spin occurs. 

In the parallel configuration and for fully polarized 
electrodes (P L = P R = 1), no reversal of the SMM's spin 
occurs and a simple linear current- voltage characteristics 
is observed. On the other hand, the linear characteristics 
disappears in the antiparallel configuration, and the cur- 
rent does not flow through the system except for the volt- 
age range where the magnetic switching of the molecule 
takes place, Fig. 3(f). 

The probabilities P m depend on the velocity c of the 
voltage increase, Eq. (7). In Fig. 4 we show (S z ) and cur- 
rent I in the parallel configuration and for several values 
of c. The magnetic switching becomes clearly visible as 
a dip in the current for larger values of c, Fig. 4(a). At 
smaller values of c, the reversal is not resolved in the cur- 
rent, Fig. 4(d). In fact, the change in c does not affect 
the time range within which the magnetic switching takes 
place, but it only modifies the dependence between the 
time and voltage scales. As a result, the transition times 
l/7m < ' + ' ) become effectively longer within the time scale 
set by the rate at which the voltage is increased. There- 
fore for the higher speeds one can observe the broadening 
of the dip. 

In summary, we showed that spin of a SMM can be re- 
versed by a spin polarized current, and the switching pro- 
cess may be visible in current when voltage is increased in 
time. Full reversal of the molecule's spin can be reached 
when at least one electrode is spin polarized. The nu- 
merical results presented above apply to the case with 
one electrode being fully spin polarized. However, the 
current-induced switching also takes place when spin po- 
larization of this electrode is smaller. The switching time 
becomes then appropriately longer. Moreover, for the pa- 
rameters assumed in numerical calculations the switch- 
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ing for positive current was only from the state | — S) 
to IS*). If the molecule would be initially in the state 
IS*), switching to the state | — S) could be achieved by 
negative (reversed) current. 
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